Orbitally controlled Kondo effect of Co ad-atoms on graphene 
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Based on ab-initio calculations we identify possible scenarios for the Kondo effect due to Co ad- 
atoms on graphene. General symmetry arguments show that for magnetic atoms in high-symmetry 
positions, the Kondo effect in graphene is controlled not only by the spin but also by the orbital 
degree of freedom. For a Co atom absorbed on top of a carbon atom, the Kondo effect is quenched 
by spin-orbit coupling below an energy scale of ~15K. For Co with spin S = 1/2 located in the 
center of a hexagon, an SU(4) Kondo model describes the entanglement of orbital moment and spin 
at higher energies, while below ~ 60 meV spin-orbit coupling leads to a more conventional SU(2) 
Kondo effect. The interplay of the orbital Co physics and the peculiar band-structure of graphene 
is directly accessible in Fourier transform tunneling spectroscopy or in the gate-voltage dependence 
of the Kondo temperature displaying a very strong, characteristic particle-hole asymmetry. 

PACS numbers: 



I. INTRODUCTION 

Graphene differs from usual metals or semiconductors 
in three important aspects: It is a truly two-dimensional 
material^ with the charge carriers resembling massless 
Dirac fermions 2 - — and the chemical potential being tun- 
able by gate voltages^. Recently, scanning tunneling 
spectroscopy experiments of graphene opened the ex- 
citing possibility to address its electronic properties lo- 
cally and to study the interaction of graphene with mag- 
netic ad-atoms. For single Co atoms adsorbed on heav- 
ily doped graphene, the observation of Kondo resonances 
with Kondo temperatures of the order of Tk = 15K has 
been reported^. Theoretically, Kondo physics in "Dirac 
materials" defined as class of materials with low energy 
Dirac type excitations has been first addressed in the con- 
text of high Tc superconductors^—. It has been demon- 
strated that even in the undoped case a Kondo effect can 
exist above a certain critical coupling between the im- 
purity spin and the Dirac electrons^. The dependence 
of the critical coupling, Kondo temperatures and impu- 
rity spectral functions on doping and localized impurity 
states has been studied in the context of graphene only 
in terms of theoretical model systems like single orbital 
Anderson models or SU(2) Kondo models^—. 

The importance of orbital physics for the Kondo ef- 
fect arises because localized spins in magnetic ions oc- 
cur almost exclusively in partially filled d or / shells. 
For graphene the same two-dimensional representations 
of the hexagonal Cq v symmetry group, that determines 
the orbital degeneracies of ad-atoms in high-symmetry 
locations, is also responsible for the band degeneracies 
in graphene at the two Dirac points. Accordingly, the 
spin of an ad-atom in the center of a carbon hexagon 
can only couple efficiently by superexchange to graphene 
bands close to the Dirac point, if it is localized in or- 
bitally degenerate levels. Therefore, the orbital degree of 



freedom and also spin-orbit coupling naturally govern the 
Kondo physics in graphene. Indeed, recent studie d 16 ' 18 ' 19 
showed that the tunneling into s-wave symmetric impu- 
rity orbitals can be strongly suppressed by graphene's 
particular symmetries, but the decisive role of the or- 
bital degree of freedom has to our knowledge not been 
studied. In general one can expect that correlation ef- 
fects will entangle fluctuating orbital and spin degrees 
of freedoms. This can lead to an SU(4) Kondo effec t 21 ' 22 . 
We show that symmetry and orbital selection rules gov- 
ern not only the coupling to the graphene bands close to 
the Dirac points but also to high-energy van-Hove sin- 
gularities. We find that virtual high-energy fluctuations 
control the size of the Kondo temperature, and, in turn, 
can lead to a strongly asymmetric gate-voltage depen- 
dence of Tk, that would be characteristic for a specific 
set of orbitals. 

The question, which model is appropriate to describe a 
certain realistic magnetic impurity system, is indispens- 
able for understanding experiments as in Ref.i but is in 
general not easy: For the classical example of Fe in Au 
studied since the 1930s, an answer could be found only 
recently 2 ^. In this article, we consider the experimentally 
important case of Co on graphene and develop a first- 
principles based model describing the Kondo physics in 
this system. This example shows, that the non-trivial or- 
bital structure of the impurity indeed controls the Kondo 
physics. 



II. DENSITY FUNCTIONAL SIMULATION OF 
CO ON GRAPHENE 

For an ab-initio description of Co on graphene we per- 
formed density functional (DFT) calculations on 6 x 6 
and 4x4 graphene supcrcclls containing one Co ad-atom 
using the Vienna Ab Initio Simulation Package (VASP)2^ 
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with the projector augmented wave (PAW ) 25 ' 26 basis 
sets. To judge the role of on-site Coulomb interaction, we 
employed a generalized gradient approximation (GGA}2! 
as well as GGA+U with U = 2eV, J = 0.9 eV and 
U = 4cV, J = 0.9 eV. We obtained fully relaxed struc- 
tures for all of these functionals. 

In agreement with Refs. [HSUli], our GGA calcula- 
tions find Co positioned above the middle of a hexagon 
on graphene (h-site) , with Co on top of carbon (t-site) or 
above a bridge site (b-site) being both more than 0.5 eV 
higher in energy. GGA predicts the electronic configura- 
tion of Co close to spin S = 1/2 at all adsorption sites. 
For U = 2eV, J = 0.9 eV h-sitc adsorption is still the 
global total energy minimum with two different electronic 
configurations of Co: First, Co can be in a spin S = 1/2 
state. The corresponding GGA+U local density of states 
at the Co site is depicted in Fig. [T^: the Co 4s-orbital 
is unoccupied in this configuration and one hole resides 
in the Co d-orbitals with El symmetry (d xz , d yz ). The 
other solution at U = 2eV, J = 0.9 cV (Fig. QJ)) yields 
the Co d-electrons can earring approximately spin S = 1 
and a state derived from the Co s-orbital directly at the 
Fermi level. We note that this solution becomes unstable 
upon decreasing of U and energetically unfavorable upon 
increasing of U. At U = 2 eV, there exists a metastable 
configuration for Co at a t-site with S = 3/2 (Fig. QJi), 
which is 0.04 eV higher in energy than the S = 1/2 h-site 
configuration. For {7 = 4 cV, J = 0.9 cV the global min- 
imum energy is found for Co with S = 3/2 on a t-site, 
which is 0.2 eV and 0.08 eV lower in energy than the h- 
and b-sites, respectively. 

In the following we will consider all relevant cases, first 
S = 3/2 on the t-site, then S — 1 on the h-site and, 
finally, S = 1/2 on the h-site as the most interesting 
case. 



III. SCENARIOS FOR THE KONDO EFFECT 

For Co adsorbed on a t-site the crystal fields have C^ v 
symmetry: the 5 d-orbitals split into two orbital dou- 
blets corresponding to the two-dimensional representa- 
tion (E) of C3 V and a singlet of the one-dimensional rep- 
resentation (Al). As can be seen from the LDOS (Fig. 
Q]b), the spin 3/2 of the Co atom is made up by a spin 
1 residing in the Co 3d-orbitals of E symmetry ferro- 
magnetically coupled to a spin 1/2 mainly from the Co 
4s orbital. The four low-energy graphene bands close 
to the Dirac points can also be decomposed into one 
two-dimensional and two one-dimensional representation 
which hybridize efficiently with the E and Al orbitals, 
respectively. In the absence of spin-orbit coupling (SOC) 
one can expect a two-stage Kondo effect: first, the direct 
coupling of the s-orbital to the C atom beneath quenches 
1/2 of S = 3/2 resulting in a spin 1 coupling to the 
two bands of E symmetry via the next-nearest C atoms 
which screen the remaining spin in a second stage. To 
estimate the strength of SOC, we calculate the magnetic 
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FIG. 1: (Color online) Orbitally resolved spin-polarized local 
density of states (LDOS) (left) and corresponding energy level 
diagrams (right) for a Co ad-atom at an h-site (a,b) and a t- 
site (c) for U = 2eV and J = 0.9 eV. At the h-site, Co has 
either spin 1/2 due to one hole in the 3d orbitals with El 
symmetry (d xz , d yz ) (a) or spin 1 with two holes in the El 
orbitals (b). For Co on a t-site we find spin 3/2 due to two 
holes in the Co 3d orbitals and one electron in the Co 4s 
orbital. 



anisotropy, -Emae = E\\ — E± ps 1.0 me V, as the energy 
difference between magnetization parallel and perpendic- 
ular to the graphene plane in GGA+U with U = 4.0 eV 
and J = 0.9 eV. This translates into different energies, 
E\s,\=i/2 — E\Sz\=3/2 ~ 1.3 mcV, depending on the z- 
component, S z , of the Co spin S = 3/2 in this config- 
uration. Ignoring Kondo physics, the ground state has 
IS, | = 3/2 and exhibits Kramers degeneracy but spin 
flips S z = 3/2 — > S z = — 3/2 induced by electron scatter- 
ing are only possible in higher order processes. Hence, 
the Kondo effect is efficiently suppressed for Co in this 
configuration, as soon as Tk 1.3meV~ 15 K. For a 
(first stage) Kondo temperature of the order of 15 K or 
above, a definite determination of the relevant low-energy 
model is very difficult, but a possible scenario is that first 
the Kondo effect partially screens 1/2 out of S = 3/2, 
then spin-orbit coupling stabilizes a low energy doublet 
(S z = ±1), which is finally screened at very low temper- 
atures by an anisotropic Kondo effect. 

For Co at an h-sitc, the situation is more interest- 
ing and a quantitative analysis easier: crystal fields of 
the Cq v symmetry decompose the d-orbitals into two 
non-equivalent two-dimensional irreducible representa- 
tions El (d xz , d yz ) and E2 (d x 2_ y 2, d xy ) plus one one- 
dimensional representation Al (d 3z 2_ r 2). For both, the 
S = 1/2 and the S = 1 configuration, the spin mainly 
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resides in the El orbitals, as shown in Fig. Q] a) and 
b). We calculated the strength of the SOC for Co on 
graphenc using VASP and obtained A = 60 meV for the 
single particle SOC operator -ffsoc ~ Al • s, with 1 and 
s being the orbital and spin angular momentum opera- 
tors, respectively. Moreover, in the GGA calculations we 
obtained the crystal field splitting from the d-level ener- 
gies as ee2 — £ei = — 0.8cV and eai — f-Ei = — 0.56eV 
for Co at the h-site. Diagonalizing a Co atom with 
S = 1 in d? configuration in this crystal field yields a 
singlet as the ground state, which is separated by about 
0.008 eV« 90 K from a doublet of first excited states. 
Hence, for a high-spin Co at an h-position an S = 1 
Kondo effect is quenched if Tk < 90 K and a much lower 
Tfc^ is only consistent with the low-spin configuration. 

For Co at an h-site with S = 1/2 in d 9 configuration 
(see Fig. [T] a) one obtains a four-fold degenerate state. 
SOC lifts this degeneracy, resulting in a twice degenerate 
atomic ground state, which is separated from a doublet of 
excited states by an energy of the order of A. In this more 
than half-filled regime, the spin- and orbital moment are 
aligned in parallel (c.f. Hund's 3rd rule). The d-holc 
resides in the highest crystal field orbitals, El, which 
have \l z \ = 1. Hence, the Zeeman splitting for out-of- 
plane magnetic fields, B z , is AE = fiBB z (gis z +gil z )/h = 
±Hb'2B z resulting in the effective g-factor of g\s z \/H = 
2. The SOC induced lifting from four fold to two fold 
degeneracy will lead to SU(4) Kondo physics above the 
scale of A and SU(2) Kondo physics at lower energies. 



IV. KONDO EFFECT OF CO AT AN H-SITE 

To address the Kondo effect in this configura- 
tion, we describe the Co at an h-site in terms of 
an Anderson impurity model: The conduction elec- 
trons residing in graphene's 7r-bands are modeled 
by a tight-binding Hamiltonian with t = — 2.97cV, 
t' = -0.073 eV, and t" = -0.33 cV quantifying the 
nearest, next-nearest and next-to-next nearest neigh- 
bor hopping, respectively^. For the Co atom, we 
consider its 3d orbitals, H imp = Y, m ,o- e \m\n m ,<r + 

^12(m,(j)^(m',(T') n m,cr n m',a' with Ti ma = d ma d m ^, 

where m is the quantum number of the z-component of 
the orbital momentum, d mj(T are Fermi operators, U is 
the local Coulomb repulsion and ei m i are the bare on-site 
energies. Here, we include hopping from the localized d- 
orbital to the nearest-neighbor C-atoms and use the Cq v 
symmetry to write the coupling of Co to graphene in the 
form 

V = J2 V \™\4n,adm,a+K.C., (1) 
m,a 

where c„ hlJ = X)<j> e% Sc j,<r/V6, Cj !<T is the Fermi oper- 
ator of electrons at carbon atom at site j and 0j is the an- 
gle between a fixed crystalline axis and the bond from site 
j to the Co impurity. All local physics, is contained in the 



local Hamiltonian, H lmp , and the hybridization function 
A mm ,(uj) defined as A mm ,(iu) = V\ m \G° mm ,(iu)V\ ml \, 
where 

G mm ,{iu) = J dfe(m|k>(m, - H^Wm') (2) 

is the bare graphene electron Green function of the states 
Cm. A mm >(iuj) = A m (iu)8 mm ' is diagonal and A m {ius) = 
A_ m (iui) by symmetry. 

The hybridization functions for different values of 
|m| = 0,1,2 are subject to selection rules imposed by 
the matrix elements (m|k): The cigenstatcs of Hq close 
to the Dirac points, K and K' = —K, transform ac- 
cording to El and E2 under Cq v with the El and E2 
being degenerate at the Dirac point. Hence, hybridiza- 
tion with to = states is cubically suppressed and 

Im A H=1 ( W )/V? = Im A M=2 H/l/ 2 2 a -7r 2ff ff 2 "J <)a 
to leading order in w. 

In contrast to the particle hole symmetry for u> — > 0, 
the hybridization functions are largely asymmetric at 
higher energies. This is caused by the El and the E2 
impurity orbitals coupling each to only one of the van- 
Hove singularities resulting from the graphene bands at 
the Brillouin zone M point: The El impurity orbitals 
as well as the graphene valence electron wave functions 
at the M point are odd under 180° rotation about the 
h-site, whereas the E2 orbitals and conduction electron 
wave functions at the M point are even under this trans- 
formation. Hence, the El hybridization exhibits a log- 
arithmic singularity, Ai(lu) ~ In |<J — Em— |j at Em- = 
t + t' — 3t" « — 2.1cV. However, there is no singularity 
in the El hybridization at the energy of the conduction 
band van- Hove singularity, Em+ = —t+t'+3t" » 1.9eV. 
For the E2 orbitals, the situation is reversed: A2(o>) ~ 
In | uj — E M + \ for w — s- E M +- 

To obtain realistic hybridization strengths, V\ and V2, 
we calculate A m by means of DFT as described in Ref. 
[HI and fit the tight-binding hybridizations via V\ and V2 
(see Fig. [5]). The tight-binding (TB) hybridizations are 
obtained in two ways: (1) by directly evaluating Eq. ©, 
which models one Co ad-atom on an infinite graphene 
sheet (TB-inf), and (2) by employing the same 6x6 
supcrccll as in the DFT calculations and performing the 
same supercell Brillouin zone integration (TB-cell). 

The high energy particle-hole / E1-E2 asymmetry is 
striking the DFT as well as in both TB hybridization 
functions. The DFT hybridization functions display 
small wiggles and the van-Hove singularities appear to 
be smeared out. Comparison of the DFT hybridization 
to the TB supcrccll hybridization shows that these two ef- 
fects are supercell artifacts. The tight-binding curves can 
be well fitted to DFT— with V\ = 1.4 eV and Fa = 1.5 eV. 
For energies above 3.4 eV and below — 3.1eV also further 
bands contribute to A m (u). They contribute to screen- 
ing and lead to a (finite) renormalization of the exchange 
coupling which we absorb in a redefinition of the bare 
exchange coupling Jo used below. In the following, we 
employ the tight-binding hybridization function for one 
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FIG. 2: (color online) Imaginary part, Im A| m |(_B), of the 
hybridization functions of the El (|m| = 1) and E2 orbitals 
(\m\ — 2) of a Co ad-atom adsorbed to a graphene h-site. 
Hybridization functions obtained from DFT and tight-binding 
(TB) models of one Co on an infinite graphene sheet (TB-inf) 
as well as the same supercell (TB-cell) as used in DFT with 
Vi — 1.4 eV and V2 = 1.5 eV are shown. 



the trend of how the Kondo temperature depends on the 
chemical potential in graphene is robust w.r.t. changes in 
Jo: Leaving Jo as fitting parameter, we predict the gate 
voltage-dependence of Tk as shown in Fig. fright. Note 
that the used values for J « leV are of order V 2 /U for 
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FIG. 3: (color online) Left: Renormalized effective coupling 
strength, J(D) pO-p)+p(m+o) = J(D)p , as function of the 
high energy cut-off D for the chemical potential [i = 0.20 eV 
and three different bare couplings Jo = 1.05, 1.00, and 
0.95 eV. Right: Kondo temperatures Tk as function of \i for 
Jo = 1.02, 1.00, and 0.98 eV. The Kondo temperatures are re- 
markably asymmetric for electron and hole doped graphene. 



Co on an infinite graphene sheet, which is free from su- 
percell artifacts. 

To estimate Kondo temperatures and their gate volt- 
age dependence, we solve the scaling equatio n 31 ' 32 

dJ(D) , . 9 , . p(u - D) + p(u + D) 

' = -N(D)J 2 (D)^ iJ-f^ '-, (3) 
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where J(D) is the renormalized exchange coupling, D the 
high energy cut-off, p the chemical potential in graphene 
and p(u>) = — Im Ai(w)/(7rl/ L 2 ). The degeneracy factor, 
N(D) = 4 for D > A and N(D) = 2 else 3 ^, accounts for 
locking the orbital- to the spin-degree of freedom below 
the energy scale of the spin-orbit coupling. 

Like all one-loop rcnormalization group equations, Eq. 
([3| is valid as long as the renormalized coupling is small, 
see Fig. [3J left. It is used to detect the energy scale 
where the strong coupling regime is approached which 
we identify with the Kondo temperature. This procedure 
correctly identifies the exponentially strong sensitivity of 
the Kondo scale on system parameters. While an exact 
calculation of the prcfactor of Tk in the limit of small 
Jo and /i ^ requires at least a two-loop calculation 3 -^, 
the one-loop equation (J3j> captures the main effect of a 
frequency-dependent density of states. Note, however, 
that the perturbative renormalization group calculation 
cannot describe the quantum-critical point£~— obtained 
for vanishing density of states as the renormalized cou- 
pling is not small in this case. 

Varying Jo for /1 = 0.2 eV (as in the experiment re- 
ported in Ref£), we find the Kondo temperature chang- 
ing by an order of magnitude for varying Jo within a few 
percent. (See Fig. [3] left.) While this hinders predic- 
tions of the absolute value of the Kondo temperature, 



realistic values of V\ and U . The remarkable asymmetry 
of the hybridization function leads to a highly asymmet- 
ric dependence of the Kondo temperatures on the gate 
voltages. If the E2 orbitals were carrying the magnetic 
moment instead of El this asymmetry would be reversed. 

Interestingly, for J > J c ~ 1.1 eV we find that the 
Kondo effect persists even for vanishing doping. This 
implies that by relative small changes (e.g. using different 
substrates) it may be possible to realize the quantum 
critical point of the pscudogap Kondo problem^— . 



V. FOURIER TRANSFORMED STM 

The symmetry of Co orbital carrying the magnetic mo- 
ment can be probed by FT-STS. In the simplest model 
(see e.g. Ref. [H), FT-STS measures the Fourier trans- 
form of the local density of states, \pk(E)\, in the vicinity 
of an impurity with the constant background of a clean 
sample being subtracted: 

Pk {E) = ~ J d 2 re lkr lm (G(r,r,E)-G°(r,r,E)). (4) 

Here, G(r, r, E) denotes the full Green function of the 
graphene-impurity system in position space represen- 
tation and G°(r,r,E) is the Green function of clean 
graphene. Using a resonant level model for the Kondo 
peak, we employ the T-matrix formalism (see e.g. Refs. 
llTll33ll34| ) in the discrete position space representation. 
Then, Eq. g]) leads to 
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FIG. 4: (color online) Fourier transformed LDOS for res- 
onances of El (left) and E2 symmetry (right) at energy 
E — 0.4 eV for a phase tv/2 of the t-matrix at this energy. 
The reciprocal lattice vectors are (2/y3, 0) and (l/\/3, 1) in 
units of 2iv/a, where a = 2.465A is the lattice constant. The 
bright spots close to the panel boundaries are centers of higher 
Brillouin zones. The insets show close-up views of the K/K' 
points which are at symmetry equivalent positions of (0, 2/3). 



where the index j labels the two atoms per graphene 
unit cell, Tj their position w.r.t. the unit cell origin and 
the k' integral extends over the first Brillouin zone. The 
Green functions occurring in Eq. ((5j) arc 2x2 matrices 
in sublattice space and obtained from the unperturbed 
graphene Green functions G k (E) by using the T-matrix: 
SG k ,, k , +k (E) = G° k ,(E)T k ,, k , +k (E)G° k , +k (E). 

In a resonant level model for the Kondo peak, we con- 
sider orbitals of El and E2 symmetry to derive the FT- 
STS patterns from Eq. using a corresponding T- 
matrix with phase 7r/2. 

The resulting Fourier transformed LDOS images are 
shown in Fig. 01 As the Kondo impurity on the h-site 
couples equally strong to both K and K' points, the inter- 
valley scattering is very strong. Due to the two sublat- 
tices, it depends however strongly on the energy E and 
the phase shift S to which extent this K-K' scattering 
leads to FT-STS intensity at the K and K' points. As 
Fig. [4] shows, there is a double arc structure of inten- 
sity around K / K' for E = 0.4 eV with the radius given 
by twice the Fermi wave vector. These structures dis- 



appear for E — > 0. The orbital symmetries manifest in 
distinct FT-STS maps with characteristic gate voltage 
and tunneling bias dependence: Upon energy, E — > —E : 
and phase shift reversal, ir/2 + -q — > 7r/2 — rj, the FT- 
STS patterns of resonances with E2 and El symmetry 
interchange. 



VI. CONCLUSIONS 

We showed that the Kondo effect of Co ad-atoms on 
graphene is controlled by the particular symmetries of 
the Co 3d orbitals originating from graphene crystal field 
splitting. Based on first-principles calculations we found 
different possible scenarios with t— or /i-sitc adsorption 
of Co and consequences for the Kondo physics. For Co 
at an /i-site we found a surprising asymmetry of Kondo 
temperatures w.r.t. the chemical potential and predicted 
characteristic FT-STS patterns. Both of these effects can 
be probed by STM. 

The importance of the orbital degree of freedom for 
the Kondo effect in graphene can be traced back to the 
symmetries underlying the peculiar band degeneracies of 
graphene at the Dirac point. Therefore the orbital degree 
of freedom is expected to control the Kondo physics in 
graphene also for other magnetic impurities occupying 
high-symmetry positions. 
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